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PREFACE 
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The  Defense  Nuclear  Agency  has  collected  and  printed  the  attached  papers 
from  the  February  25-27  1986  Global  Effects  review  as  a  service  to  the 
community.  The  Defense  Nuclear  Agency  takes  this  opportunity  to  express  Its 
gratitude  to  the  numerous  participants  In  the  Global  Effects  review. 

The  technical  papers  enclosed  Include  all  those  which  were  received  by  DNA 
prior  to  the  closing  date  of  28  April  1986.  Where  papers  are  missing  their 
place  Is  occupied  by  the  abstract  received  prior  to  the  meeting. 

The  Inclusion  of  a  paper  In  this  proceeding  does  not  necessarily  imply 
endorsement  of  the  results  of  the  research  reported  or  conclusions  which  might 
be  drawn  from  that  research.  It  is  the  opinion  of  the  Defense  Nuclear  Agency 
that,  while  good  progress  is  being  made  in  Improving  our  understanding  of 
Global  Effects,  the  results  to  date  are  tentative  and  preliminary  and  should 
not  be  used  for  planning  beyond  the  planning  of  future  research. 
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LABORATORY  STUDY  OF  WET  COAGULATION 


GLOBAL  EFFECTS  PROGRAM  TECHNICAL  MEETING 


FEBRUARY  26,  1986 


GEORGE  CARRIER* 
FRANK  FENDELL 
DAN  KWOH 
BRUCE  LAKE 


•Harvard  University,  Cambridge,  Massachusetts  02138 


TRW  SPACE  AND  TECHNOLOGY  GROUP 
ONE  SPACE  PARK 

REDONDO  BEACH,  CALIFORNIA  90278 


OBJECTIVE 


Do  laboratory  experiment: 

(1)  To  see  if  light  extinction  by  soot 
particles  changes  after  confinement  in 
a  cloud 

(2)  To  measure  size  distribution  and  number 
density  of  soot  particles  before  and 
after  confinement  in  a  cloud 


r 

APPROACH 


(1)  Introduce  soot  particles  and  moist  air  into 
cloud  chamber 

(2)  Measure  extinction,  number  density  and  size 
distribution  of  soot  particles 

(3)  Expand  cloud  chamber  adiabetical ly  to  form  cloud 

(4)  Sit  for  sometime  for  scavenging,  wet  coagulation 
to  take  place 

(5)  Compress  cloud  chamber  to  evaporate  cloud 

(6)  Measure  extinction,  number  density  and  size 
distribution  of  soot  particles  again 
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OUTER  WALL 


INNER  WALL  -  AI 
SEAL 


REFRIGERANT 


PARTICLE 

EXTRACTION 

PROBE 


OPTICAL  WINDOW 


CLOUD  CHAMBER:  APPROACH 


•  Wall  cooling  by  simply  passing  refrigerant  through 
at  controlled  cold  temperature  and  flow  rate.  Cool 
from  R.T.  to  ~  -60°C  in  1/2  hour. 

•  Gas  cooling  by  adiabetic  expansion  using  piston. 
Expansion  ratio  ~  1:3 

•  Minimize  leakage  by  piston  by  evacuating  and  equalizing 
pressure  on  back  side  of  piston 

•  Monitor  wall  and  gas  temperature  by  thermocouple  or 
transistor  thermometer 

•  Make  gas  temperature  to  follow  wall's 

•  Aim  for  wall  0.1-0.3eC  warmer  than  gas  to  prevent 
condensation 

•  Make  top  of  chamber  warmer  than  bottom  to  suppress 
convection 


EXTINCTION  MEASUREMENT:  SCHEMATICS 


EXTINCTION  MEASUREMENT:  APPROACH 


.  lO^-lO5  0.1  fi  particles/c.c,*=H>  10'^-10"3/m  extinction 
coefficient 

•  Use  "laser-cavity  extinction"  method: 

-  make  soot-laser  air  part  of  laser  cavity 

-  laser  output  highly  non-linear  function  of 
cavity  loss 

-  equivalent  measurement  pathlength  ~  kilometers 

-  use  variable  Brewster  angle  flat  as  in-cavity 
calibration  device 

-  evacuate  chamber  either  before  or  after  dust 
measurement  for  calibration 
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LASER  OUTPUT  VS.  CAVITY  LOSS 
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FIGURE  2.  Visual  range  and  output  power  versus  percent  cavity  loss.  The 
prototype  uses  a  IOO-cm  sample  cavity  length.  With  current  state-of-the-art 
high  reflectivity  mirrors,  the  sample  cavity  length  can  be  reduced  to  less  than 
30  cm  with  a  Q  of  2,000  to  3,OOQ. 
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BREWSTER  ANGLE  FLAT  AS  IN-CAVITY  CALIBRATION  STANDARD 


Incident  Angle 

LOSS 

56.0 

4  x  10'7 

56.5 

6.4  x  10'5 

57.0 

2.4  x  10'4 

57.5 

5.3  X  10"4 

58.0 

9.7  X  }0'4 

58.5 

1.53  x  10~3 

59.0 

2.25  x  10'3 

60.0 

4.18  x  10~3 
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SOOT  PARTICLE  NUMBER  AND  SIZE  MEASUREMENTS 


EXTRACTION  AND  ELECTRON  MICROSCOPY 


r 


Methods  of  Extraction: 

(1)  Impaction  — many  lose  small  particles  X 

(2)  Suctioning  and  filtration  —  agglomeration  may 

occur  in  suctioning  probe  or  filter?  x 

(3)  "Grabbing"  and  gravity  settling  —  disturb 
particles  the  least,  0.1m  particles  settle 

in  3  hours  V 

Methods  of  Electron  Microscopy: 

(1)  Scanning  electron  microscopy  (SEM)  —  0.1  m 
resolution  may  be  difficult,  need  to  coat 
particles  with  100-300  A  gold  layer? 

(2)  Transmission  electron  microscopy  (TEM)  —  sample 
preparation  more  involved  but  well  developed. 

Mount  Collidian  film  on  screen  or  slide. 


OPTICAL  METHODS  AVAILABLE  FOR  MEASURING  NUMBER  DENSITY 
AND  SIZE 


Ensemble  Measurement 
(beam  measurement) 


Scattering  extinction  ratio 
2  angle  intensity  ratio 
diffusion  broadening  spectroscopy 
2  wavelength  extinction  ratio 
polarization  ratio 


Particle  Measurement 
(focused  measurement) 


2  angle  intensity  ratio 

2  beam  interferometry 

white  light  integrated  forward 
scattering 


COMMON  SOURCES  OF  ERROR  FOR  ALL  OPTICAL  TECHNIQUES 


1.  Non-spher icitv 

maximum  error  -  A0% 

2.  Index  of  Refraction  Uncertainty 

for  highly  absorptive  particles  like  sooty 
a  few  percent 

3.  Depth  of  Field  (for  focused  measurements) 

broadens  size  distribution,  <  30% 


ft 


OPTICAL  METHODS  SELECTED 


For  small  particles  (<  0,1  n)  —  extinction/scattering 

ratio  or  2  angle  intensity 
ratio  (large  angle  pair) 

For  large  particles  (>  several  ^ )  —  2  angle  intensity 

ratio  (small  angle  pair) 
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ME  SCATTERING 


SUMMARY  OF  WHAT  WE  PROPOSE  TO  DO 


(1)  Measure  extinction  accurately 

(2)  Measure  particle  size  distribution  and  number 
density  accurately  by  extraction  and  TEM 


(3) 


Use  other  optical  techniques  for  real  time 
monitoring  and  as  secondary  measurements 
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MICROPHYSICAL  PROPERTIES  AND  REMOVAL  PROCESSES  OF  SOOT 


John  Hallett,  Barry  Gardiner,  Dennis  Lamb, 
Fred  Rogers,  Richard  Pitter,  Jim  Hudson 


Desert  Research  Institute 
Reno,  Nevada  89506 


Viewgraphs  -  February  1986 
NASA  AMES  DNA  Meeting 


NUMBER  OF  DROPS  FROZEN  PER  DEG  C 


TEMPERATURE  (DEG  C) 


Figure  10.  Number  of  drops  frozen  per  degree  temperature 
change  as  a  function  of  temperature. 
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FIG.  12:  COAGULATION  RATE  OF  AEROSOL  WITH  CLOUD 
DROPLETS,  TIME  CONSTANT  ^lOOOS 
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:  IGURE  17: 


Figure  18  a)  Cloud  produced  on  natural  aerosol  In  an  expansion  cloud 
chamber. 
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b)  Cloud  produced  on  acetylene  soot  in  an  expansion  cloud 
chamber. 
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SOOT  PARTICLE  DISTRIBUTION  C23-25  OCT  ' 85) 
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t  =  200  s 


CRITICAL  RADIUS 


CRITICAL  SUPERSATURATION,  % 

FIG.  24:  MODEL  DROPLET  GROWTH  10®  cm'3 
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FIG.  25:  MODEL  DROPLET  GROWTH  100  cm-3 


—  SOLUBLE 

-  INSOLUBLE 


TOTAL  CONCENTRATION,  cm 


FIG.  26:  MODEL  NUMBER  AND  FRACTION  ACTIVATED 


CONCLUSIONS  FROM  OUR  STUDY: 

1) .  Lab/field  soot  complex;  0.1  to  100  pm. 

2) .  Lab  soot  acetylene/crude  oil  CCN  ratio  few  to  10%  increasing  with  time 

CN 

3) .  Ice  nucleation  unexciting. 

4) .  Soot  collects  on  outside  of  drops  mm  -  50  um  as  evaporation  proceeds. 

(i)  Independent  particles 

( i i )  lily  pads 

(iii)  solidified  surface 

(iv)  collapsed  soccer  ball 

(v)  surface  resuspension 

5) .  Sooty  cloud  doesn't  change  extinction  over  lO^s  by  scavenging. 

6) .  CCN  -  cloud  droplets  -  ice  removal  or  coalescence;  only  partial 

pluvial  constipation. 


FUTURE  WORK: 

i.  evaporated  cloud  droplets  and  soot  coagulation. 

ii.  scattering  +  extinction  of  sooty  cloud. 

iii.  role  of  qrowth/evaporation  in  scavenging  process  (a)  droplets, 
(b)  ice  * 

iv.  soot  in  a  supercooled  cloud  -  rime  growth. 

FIELD: 

i.  soot  from  oil  fire.  CN,  CCN,  replica. 

ii.  cloud  droplets  on  replica,  modification  and  scavenging  in  cap 
cloud.  (Lodi?) 

iii.  deeper  cloud  and  ice  role  in  scavenging.  Ice  from  repl ica 
(Chapleau? ) 


MAJOR  PROBLEM: 

Firestorm  dynamics  to  give  right  microphyics  for  maximum  scavenging. 


Chemical  Scavenging  in  the  Atmosphere:  Smoke-Ozone  Reactions- 


N.  deHaas,  R.  M.  Fristrom,  M.  J.  Linevsky  and  D.  M.  Silver 
Applied  Physics  Laboratory,  The  Johns  Hopkins  University 
Laurel,  Maryland  20707 


ABSTRACT 

A  laboratory  apparatus  has  been  constructed  to  examine  the  reaction 
of  smoke  particles  with  ozone.  A  stopped-flow  experimental  procedure 
has  been  developed  that  permits  (1)  a  measurement  of  the  total  extinc¬ 
tion  of  a  He-Ne  laser  beam  by  ozone ( Chappuis  band)  and/or  by  various 
sub-micron  smoke  particles;  (2)  a  measurement  of  the  right-angle  scat¬ 
tering  of  a  Nd-Yag  laser  by  the  smoke  particles;  and  (3)  a  sampling  of 
particulates  for  examination  under  the  scanning  electron  microscope. 
Ozone  (without  smoke)  is  stable  in  the  apparatus  over  a  period  of  hours; 
smoke  (without  ozone)  is  also  stable  over  a  period  of  hours.  Smoke- 
ozone  mixtures  begin  decaying  immediately  and,  under  the  experimental 
conditions  employed,  they  have  half-lives  of  a  half-hour  to  two  hours 
depending  on  the  fuel  source  of  the  smoke  particles.  Specifically,  the 
half-life  is  defined  to  be  the  time  in  which  the  optical  extinction  at 
6328  A  is  reduced  to  one-half  its  initial  value.  Translating  these 
results  to  "nuclear  winter"  atmospheric  conditions  yields  half-lives  in 
the  range  of  a  month.  These  experiments  provide  solid  qualitative 
evidence  for  the  existence  of  atmospheric  chemical  scavenging  of  smoke 
and  ozone. 


This  work  was  supported  in  part  by  the  Naval  Sea  Systems  Command  under 
Contract  N00024-85-C-5301 . 


Submitted  for  presentation  at  the  DNA  Global  Effects  Program  Technical 
Meeting,  25-27  February  1986  by  D.  M.  Silver,  (301)  953-6265. 
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HETEROGENEOUS  REACT I ONS 
IIM  A  SOOT  LADEN  ATMOSPHERE 

DAV I D  M .  GOLDEN 
M  I  CHEL  Or-  ROSS  I 

LUDWIG  BROUWER 
SHERRY  STEPHENS 

DEPARTMENT  OF  CHEM I CAL  KINETICS 
CHEMICAL  PHYS I OS  LABORATORY 
SRI  I NTERNAT I ONAL 


Table  1.  Collislonal  reaction  probabilities  on  a 
H2S04  surface  at  300  K. 


Collisional 

reaction 

Species  probability 


M 


h202 

7  .8 

X 

10“4 

HNOj 

£ 

2.4 

X 

10-4 

ho2nd2 

2.7 

X 

10-s 

CIONOj 

1  .0 

X 

io-5 

n205 

£ 

3.8 

X 

IO-5 

h2o 

~ 

2.0 

X 

io-3 

Nil  3 

> 

1  .0 

X 

10_s 

OH 

4  .9 

X 

io-4 

o3 

< 

1  .0 

X 

1 0“* 

NO 

< 

1  .0 

X 

io-* 

NO, 

< 

1  .0 

X 

io-6 

9°, 

< 

1  .0 

X 

io-8 

Alkenes 

< 

1  .0 

X 

io-8 

Alkanes 

< 

1  .0 

X 

10-e 

cf4 

< 

1  .0 

X 

io-8 

cci2f2 

< 

1  .0 

X 

v 

1 

O 

r-4 

0(3P) 

< 

1  .0 

X 

IO"8 

N 
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1  .0 
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10'8 
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PROPOSED  MODELING  STUDIES  ON  REMOVAL  BY  CLOUD  PROCBSS 
OP  MASSIVE  EMISSIONS  OF  PARTICLES 


In  Young  Lee 

Atmospheric  Physics  Program 
Environmental  Research  Division 
Argonne  National  Laboratory,  Argonne,  IL  60439 

ABSTRACT 

In  general,  atmospheric  aerosols  originate  by  several 
processes:  (1)  condensation  and  sublimation  of  vapors  and  the 
formation  of  smokes  in  combustion,  (2)  reactions  among  trace 
gases,  (3)  mechanical  disruption  and  dispersal  of  matter  at  the 
earth's  surface  as  mineral  dusts,  and  (4)  coagulation  of  nuclei 
to  produce  larger  particles  of  mixed  constitution.  The  particles 
participate  actively  during  cloud  formation  either  as  cloud 
condensation  nuclei  or  ice-forming  nuclei,  or  passively  by  being 
removed  by  cloud  droplets  and  raindrops  through  processes  of 
collision  and  coalescence.  The  evolution  of  dust  and  smoke 
particles  is  additionally  controlled  by  dynamio  processes  such 
as  transport,  turbulent  mixing  and  deposition.  These  natural 
processes  of  aerosol  evolution  would  be  greatly  affected  if 
extraordinary  amounts  of  dust  and  smoke  particles  were  injected 
into  the  atmosphere  by  blasts  or  large  fires. 

The  scavenging  efficiency  for  aerosols  produced  by  blasts 
and  fires  can  be  examined  in  numerical  simulations  which 
describe  coagulational  processes  involving  aerosols,  cloud 
droplets  and  raindrops.  The  investigation  of  the  evolution 


within  clouds  of  dust  and  smoke  particles  of  very  high 
concentrations  is  important  for  clear  understanding  of  the 
resulting  global  budget  of  aerosols  and  for  better  assessment  of 
the  potential  climatic  impact.  The  collision  efficiency  between 
partioles  can  be  formulated  from  a  combination  of  theoretical 
analyses  and  available  experimental  data.  A  numerical  module 
for  wet  removal  can  be  used  to  evaluate  scavenging  efficiency, 
modification  of  cloud  properties,  and  applicability  of  cloud 
microphysics  parameterizations  currently  in  use. 

A  cloud  model  originally  developed  at  Argonne  to  study 
cloud  microphysics  parameterizations  has  been  modified  to 
examine  the  behavior  of  submicron  aerosols  in  clouds . 
Lagrangian  computations  following  a  rising  parcel  of  air  were 
made  in  order  to  understand  and  identify  mechanisms  that  control 
the  fate  of  smoke  and  dust  particles  in  clouds.  The  air  parcel 
was  assumed  to  be  saturated  initially  and  to  contain  aerosols  of 
high  number  density  (1650  cm-3).  Some  portion  of  the 

aerosols  in  the  parcel  could  be  activated  to  form  cloud  droplets 
and  the  rest  could  remain  inactivated  due  mainly  to  the  Kelvin 
effect.  Thereafter,  the  spectrum  was  allowed  to  evolve  and  form 
raindrops.  The  model  treated  particle  sizes  ranging  from  0.04 
fx m  to  164  fim  in  radius. 

Computations  have  been  made  for  40  min  of  simulated  time. 
The  spectral  evolution  and  the  role  of  each  microphysical 
mechanism  in  the  spectral  evolution  are  presented  in  Figs.  1  and 


2,  respectively.  Figure  1  shows  that  the  . peotrum  evolves 
quickly  during  the  first  five  minutes  of  cloud  development  due 
mainly  to  condensation,  which  affeots  the  growth  of  particle 
sizes  approximately  up  to  the  mode  radius  near  10  /xm.  The 
spectral  evolution  starts  to  produce  raindrops,  associated  with 
sizes  greater  than  about  50  1b  radius,  after  10  min  of 
simulated  time.  The  submicron  size  particles  are  not  removed  or 
activated  until  simulated  time  exceeds  20  min,  but  thereafter 
they  are  quickly  scavenged  as  the  cloud  starts  to  accumulate 
drops  greater  than  about  100  /x m  in  radius.  In  Fig.  2,  we 
see  clearly  that  the  combined  processes  of  condensation  and 
coagulation  can  remove  submicron  partioles  much  more  effectively 
than  either  process  alone.  We  also  found  that  the  removal  of 
submicron  particles  is  augumented  if  precipitation  is  included. 

Figure  1.  Number  distributions  of  aerosols,  cloud  droplets 
and  raindrops  after  various  periods  of  simulation. 

Figure  2.  The  role  of  cloud  microphysical  processes  of 
condensation  and  coagulation  in  the  spectral  evolution  of 
particles . 
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HE  UMK  CLOUD  SIMULATION  CHAMBERS 


The  UMR  cloud  simulation  program  is  based  on  two  cooled-wall 
expansion  chambers  designed  to  produce  uniform  cooling  comparable 
to  that  produced  in  cloud  updrafts.  The  smaller  chamber,  called 
Proto  II,  is  currently  producing  cloud,  and  the  larger — more 
sophisticated — Romulus  chamber  is  scheduled  for  completion  in  the 
summer  of  1985. 

These  chambers  feature  thermoelectric  modules  to  cool  the 
walls  at  the  same  rates  as  the  expansion  cools  the  gas.  Thermal 
wall  effects  are  thus  effectively  removed  and  a  closed  parcel 
thermodynamic  environment  similar  to  that  found  in  the  core  of  a 
convective  cloud  in  the  absence  of  turbulent  mixing  can  be 
simulated . 

After  a  sample  of  moist  aerosol-laden  air  has  been 
introduced  into  the  chamber,  the  sensitive  volume  undergoes  a 
controlled  expansion  which  cools  the  gas  and  produces  the 
increased  relative  humidity  required  for  cloud  formation.  The 
chamber's  initial  condition  can  be  set  at  any  desired  state  of 
relative  humidity  less  than  100%.  The  cooling  capability  is 
sufficient  to  allow  ice  nucleation  and  crystal  growth  studies  to 
be  carried  out. 

The  creation  of  a  viable  cloud  simulation  program  has 
necessitated  developmental  work  on  a  broad  front  so  that  the 
required  science  and  technology  would  come  together  to  support 
the  experimental  effort  at  about  the  same  time.  Much  of  this 
effort  has  gone  into  the  development  of  the  lollowing  peripheral 
and  supporting  equipment: 
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1.  Air  preconditioning  system  for  drying  and  filtering  large 
volumes  of  air  (see  Fig.  A2)  for  flushing  the  chamber  and 
for  use  by  aerosol  generation  equipment. 

2.  Precision  saturator  columns  and  thermal  regulation  columns 
for  relative  humidity  (see  Fig.  A3). 

3.  A  variety  of  aerosol  generation  techniques  have  been 
developed  for  the  production  of  artificial  nuclei.  Much 
effort  has  gone  into  characterizing  aerosols  and  rendering 
them  monodispersed.  Natural  nuclei  can  be  employed. 

4.  Precision  continuous  flow  thermal  diffusion  chamber  (CFDs) 
for  characterizing  the  critical  activation  supersaturation 
spectrum  of  the  conder nation  nuclei  have  been  developed. 
Three  vertical  plate  CFDs  are  currently  operational.  Three 
vertical  plate  CFDs  are  currently  under  construction.  These 
will  operate  simultaneously,  two  in  the  CFD  mode  and  one  in 
the  isothermal  mode,  to  cover  the  entire  CCN  spectrum  with 
high  resolution  and  will  provide  automatic  data  acuisition 
in  a  format  for  immediate  assimilation  by  the  simulation 
chamber  minicomputer  system  to  aid  in  establishing  last 
minute  parameters  for  experiment  control. 
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5.  Optical  and  electrical  aerosol  characterization  equipment 
for  measuring  the  size  distribution  of  both  wet  and  dry 
aerosol  particles.  A  GE  CNC  counter  and  a  UMR  Aitken 
nucleus  counter  are  also  available  for  mesuring  total 
particle  concentrations  and  sensing  the  Aitken  nuclei 
content  of  sample  air.  This  makes  possible  scavenging 
experiments. 
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6.  Laser  light  scattering  equipment  enables  one  to  determine 
the  Mie  scattering  intensity  at  one  angle.  The  total 

t 

attenuation  of  a  light  beam  (from  a  light  emitting  diode 

operating  at  a  different  wavelength  from  the  argon  ion 

laser)  is  also  measured  with  a  precision  of  three  parts  in 
4 

10  .  A  photograph  is  also  taken  with  flash  illumination  in 
order  to  verify  the  droplet  concentration  at  specified 
intervals.  These  optical  systems  are  described  in  more 
detail  in  Appendix  C. 

7.  The  digital  data  acquisition  and  control  system  is  built 
around  a  NOVA  840  and  a  NOVA  3  minicomputer.  When  both  are 
operational,  a  greater  real  time  computational  capacity  is 
available  for  experiments.  State-of-the-art  techniques  are 
employed  extensively  through  the  system,  particularly  in  the 
measurement  and  control  of  temperature  and  pressure. 

8.  A  wide  variety  of  simulation  chamber  control  options  have 
been  developed.  These  can  be  readily  configured  to  best 
suit  the  demands  of  individual  experiments.  It  is  possible 
to  run  various  microphysical  models  in  real  time  on  one  of 
the  minicomputers  and  use  the  output  in  the  control  of  the 
chamber  or  to  compare  with  the  output  of  the  various  optical 
sensors.  In  this  area  we  envision  a  continuous  evolution  of 
capability  as  new  experiments  are  developed. 
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UMR  CLOUD  SIMULATION  FACILITY 
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TABLE  I 

PROTO  II  CHAMBER  SPECIFICATIONS 


?emperature 


Prototype 

Chamber 


Range 

Control 

Holding:  At  fixed  temperature 

RMS  deviation  from 
spatial  average,  40  sensors 

H  U  UU  V* 

0 . 030°C 

Cooling  Mode:  Maximum  rate,  at  least 
Wall  temperature  lag 
from  control  point 

RMS  deviation  from 
spatial  average 

10°C/min. 

0.01°C/{ “C/min. )* 

0. 080°C  @  6°C/min. 
0. 040°C  @  3°C/min . 

Wall  Measurement  Precision 

+0.002°C 

Pressure 

Reference  pressure  absolute  accuracy 

+0.01%,  dead 
weight  pressure 
gauge 

Dynamic  sensor  range 

0  to  100  Torr 
differential 

Resolution 

+0.05  Torr 

Control : 

Holding:  Standard  deviation  (time) 

Expansion:  Maximum  rate 

0.5  Torr 

3.9  Torr/sec 

12  (C/min.  wet 
ad iabatic ) 

Offset  from  command  valve 

Standard  deviation  with  (time) 

0.5  Torr 

0.36  Torr 

Wall  Cooling  Method 

Thermoelectric 
modules  plus  fluid 
heat  exchange  in 
heat  sink 

Volume 

48  cm  dia  x 

61  cm  ht 

OR 

48  cm  dia  x 

122  cm  ht 

•Cooling  Rate 

$ 
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Scavenging  of  submicron  aerosols  by  clouds. 

R. F. Pueschel,  NASA  Ames  Research  Center,  Moffett  Field,  CA 
94040. 

Scavenging  by  clouds  and  precipitation  is  the  most  effective 
mechanism  by  which  aerosols  are  removed  from  the  atmosphere. 
Data  on  the  relationships  between  aerosol  and  cloud  spectra 
and  concentrations  before,  during  and  after  cloud  formation 
show  the  extent  and  time  scales  with  which  the  aerosol-cloud 
interaction  occurs.  Specifically,  the  following  results  will 
be  discussed: 

(1)  The  accumulation  mode,  the  longest  lived  portion  of 
the  atmospheric  aerosol*,  provides  the  cloud  condensation 
nuclei  around  which  drops  form.  Particles  as  small  as  0. 1  urn 
diameter  are  ingested  into  cloud  drops  during  the  cloud's 
formation;  the  nucleation  scavenging  of  submicron  particles 
is  the  more  efficient  the  larger  their  size. 

(2)  In  polluted  airmasses,  there  are  more  submicron 
particles  removed  from  the  aerosol  than  there  are  cloud 
drops  formed.  This  is  different  from  clean  air  situations 
where  a  1:1  relationship  between  the  loss  of  submicron 
particles  and  the  formation  of  cloud  drops  exists.  Thus,  in 
polluted  air  a  mechanism  in  addition  to  nucleation 
scavenging  has  to  be  invoked  to  account  for  the  aerosol 
losses.  Stephan  flow,  a  hydrodynamical  movement  of  a  vapor- 
gas  mixture  normal  to  the  surface  of  a  condensing  liquid,  to 
compensate  for  gas  diffusion  away  from  the  surface,  is  a 
possible  explanation. 

(3)  Nucleation  scavenging  during  cloud  formation  is  180 
times  faster  than  is  in-cloud  attachement  of  the  residual 
aerosol  after  the  cloud  has  formed. 
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THE  ROLE  OF  BUOYANCY  INDUCED  TURBULENCE 
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Applications  International 


Science  Applications  International  Corporation  -  - 
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Single  Plume  Characteristic  Regions 
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Strong  Convective  Mixing 
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Rume  Structure 


BUOyANCy  GENERATED  TURBULENCE 


Science  *fpHr^hemt  tntomuL  wti  Corporation 


Science  Applications  International  Corporation - 

EFFECT  OF  MODELING  ASSUMPTIONS  ON  PREDICTION  OF 
CENTERLINE  TEMPERATURE,  HOT  BUOYANT  JET 
(Uj  =  21.6  ft/sec) 
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•  Scaling 

•  Multiple  Plume  Interactions 

•  Fuel  Effects 


On  the  Scavenging  of  Aerosol 
Particles  by  Snow  Crystals 


Pao  K.  Wang 
Department  of  Meteorology 
University  of  Wisconsin 
Madison,  WI  53706 


Aerosol  Removal 
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Wet  Removal 
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Scavenging 


Below-Cloud 

Scavenging 
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Previous  Research  In  Precipitation  Scavenging 


1.  Rain  Scavenging 


Wang,  P.K.,and  H.R.  Pruppacher . 1977b 

"An  experimental  determination  of  the  efficiency  with  which  aerosol 
particles  are  collected  by  water  drops  in  subsaturated  air" 

J.  Atmos.  Sci ■ ,34, 1664-1669. 

Wang.P.K. ,S.N.  Grover,  and  HR.  Pruppacher . 1978 
"On  the  effect  of  electric  charges  on  the  scavenging  of  aerosol 
particles  by  clouds  and  snail  raindrops" 

J.  Atnos .  Scl ■ .35. 1735-1743. 

Wang,  P.K.,  1979a 
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Theoretical  Consideration 


Basically  a  convective  diffusion  problem 
The  mechanisms  involved  are: 

1.  Inertial  impaction 

2.  Brownian  diffusion 

3.  Thermophoresis 

4.  Diffusiophresis 

5.  Electric  forces 

—Static  charges 
-Extrenal  electric  field 

6.  Turbulence 


Models: 

1.  Trajectory  Model  -  Brownian  diffusion  not  treated 

valid  for  particle  radius  >  0.5  micron. 

2.  Flux  Model  -  inertial  impaction  not  treated 

valid  for  particle  radius  <=  0.5  micron. 
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PHYSICAL  SIMULATION  OF  LARGE- AREA  FIRE  PLUMES 

by 

M.  Poreh,  J.  E.  Cermak  and  J.  A.  Peterka 

ABSTRACT 

The  dynamics  of  large-area  fire  plumes  has  been  investigated 
by  physical  simulation  using  small-scale  models.  The  models  use 
heat  or  helium  to  generate  the  buoyancy  flux.  They  neglect  the 
radiative  heat  transfer  and  assume  a  relatively  shallow 
combustion  layer,  an  assumption  which  is  found  to  be  consistent 
with  the  measurements  in  the  model,  but  maintain  the  correct 
dimensionless  parameters  which  determine  the  dynamics  of  the 
plume  outside  the  combustion  layer. 

Four  regions  of  distinct  flow  characteristics  have  been 
identified.  A  combustion  layer  whose  height  is  small  compared  to 
the  horizontal  dimension  of  the  fire.  A  radial  convective 
boundary  layer,  where  radially-inward  induced  winds  carry  the 
rising  plume  toward  the  center.  A  central  core  and  a  convective 
plume  above  it,  which  rises  to  large  heights. 

Vigorous  mixing  with  the  ambient  air  occurs  in  the  radial 
convective  boundary  layer,  where  the  hot  burned  fuel  rises  in  the 
form  of  thermals.  Such  thermals  are  formed  for  both  continuous 
combustion  layers  and  for  large  fires  which  are  made  of  many 
distinct,  close,  small  fires.  It  is  found  that  these  two  cases 
produce  similar  plumes. 

The  simulations  show  that  the  plume  above  the  burning  area 
contracts  considerably  as  it  rises,  and  that  the  upper  convective 
plume  behaves  as  a  weak  plume  with  a  virtual  origin  below  the 
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surface,  at  Z/R  of  the  order  of  -1.5.  The  subsequent  plume 
rise  in  an  upper  stably-stratified  atmosphere  would  thus  be  much 
smaller  than  that  of  a  plume  originating  from  a  point  source  at 
the  ground,  with  the  same  buoyancy  flux. 

The  simulations  also  show  that  the  induced  radial  velocities 
of  a  large-area  fire  in  a  calm  neutral  atmosphere  are  not  as 
large  as  those  estimated  by  previous  numerical  and  analytical 
models.  In  addition,  the  study  suggests  that  the  kinetics  of  a 
single  fire,  which  is  surrounded  by  numerous  other  fires  would 
not  be  much  different  from  that  of  an  individual  fire  in  a  calm 
atmosphere,  as  oxygen  depletion  is  found  to  decrease  with  the 
size  of  the  fire. 
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•  The  simulation  is  based  on  the  conventional  approach  to 
fluid  modelling.  All  the  important  dimensionless  parameters 
and  boundary  conditions  must  be  the  same  in  the  model  and 
prototype  (njn  =  np) . 

•  The  Reynolds  number  in  the  model  cannot  be  matched  with  that 
of  the  prototype.  However,  it  is  shown  to  be  large,  so  that 
'Reynolds  number  independence'  is  expected. 

•  The  models  neglect  radiative  effects. 

•  It  is  assumed  that  the  relative  thickness  of  the  combustion 
layer  (CL)  in  large-area  fires  is  small  compared  to  the  size 
of  the  fire.  Thus,  the  model  simulates  the  flow  above  the 
relatively  shallow  CL. 

•  The  flow  above  the  CL  is  assumed  to  be  governed  by 

a  -  The  buoyancy  flux  per  unit  area  produced  in  the  CL, 
which  is  proportional  to  the  burning  ratio  (j/sec)  [a] 
=  l2/t3 . 

R  -  The  size  (radius)  of  the  (assumed  circular)  CL. 

-  Dimensionless  boundary  conditions  describing  the 
ambient  wind,  ambient  circulation  and  atmospheric 
stratification. 

The  presentation  focuses  on  Large  Area  Fires  (LAP)  in  a 
calm,  neutrally  stratified  atmosphere.  It  has  been  shown, 
however,  that  the  physical  simulation  can  include  the  effects  of 
ambient  stratification  and  circulation  (Poreh  et  al. ,  1986). 


190 


The  above  assumptions  imply  that  the  velocity  field  u(x)  and 
the  effective  gravitational  field  g'(x)  =  gAp(x)/pa  are  given  by 


=  f  ( x  /  R) 

OR)1'3 


(1) 


and 


_ SLlUi 

a1'3  r-1'3 


g (x/R) 


(2) 


Buoyancy  flux  in  the  models  have  been  generated  by: 

(1)  Heat  flux,  as  suggested  by  R.  R.  Long. 

(2)  Flow  of  lighter-than-air  gas  (helium)  from  a  porous 
circular  plate. 

It  is  noted  that  the  value  of  Ap/pa  Qf  helium  is  approximately 
equal  to  that  of  burned  fuel  at  stoichiometric  proportions 
(BFSP) .  Namely,  ideal  burning  of  the  fuel  without  mixing  with 
surplus  air.  Thus,  the  model  simulates  an  area  source  of 
buoyancy  produced  by  a  flux  of  BFSP  before  it  mixes  with 
additional  air. 

Theoretical  considerations  and  previous  experience  suggest 
that  Ap/ pa  of  the  buoyant  gases  at  the  source  need  not  be  the 
same  in  the  model  and  the  prototype,  unless  a  strong  whirl  is 
expected.  This  assumption  is  confirmed  by  the  experimental 
results  (Fig.  6).  The  value  of  Ap/pa  in  one  of  the  model  (R  =  50 

a 

cm)  was  matched  with  the  expected  value  in  a  R  =  2  km  urban  fire. 


•  Two  circular  helium  models  with  R  *  27  cm  and  R  *=  SO  cm  were 
used. 

•  A  two-dimensional  area  source  (2  m  x  2  m)  was  placed  in  a  2 
m  wide  test  section  of  a  wind  tunnel. 

Summary  of  the  Experimental  Results 

•  The  shape  of  the  plume  based  on  flow  visualization  and 
helium  concentration  measurements  is  shown  in  Figures  2  and 
3.  Considerable  necking  of  the  plume  is  observed,  similar 
to  the  necking  of  a  14  acre  fire  in  Project  Flambeau  (Figure 
4) . 

•  Four  regions  of  distinct  flow  were  identified  (Figure  3): 

(1)  A  relatively  shallow  CL. 

(2)  a  radial  convective  boundary  layer. 

(3)  A  central  core. 

(4)  A  convective  plume.  The  convective  plume  which  can  be 
described  as  a  plume  from  a  virtual  point  source  at 
approximately  zy  -  -1.5  R. 

The  assumption  of  a  relatively  small  thickness  of  the 
combustion  layer  is  based  on  the  flame  height  correlation  of 
Thomas  (1963)  shown  in  Figure  5.  The  value  of  m'  ' /  ( pQ  \|gD)  for 
large  urban  fires  is  expected  to  be  of  the  order  of  10-5,  so  that 
the  ratio  of  the  flame  height  to  the  size  of  the  fuel  bed  (L/D) 
is  expected  to  be.  according  to  this  correlation,  of  the  order  of 
io“2. 
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The  model  does  not  have  a  CL.  However,  a  'pseudo  combustion 
layer'  (PCL)  has  been  defined  by  the  height  where  the  injected 
helium  mixes  with  6  times  air  by  mass.  This  is  the  order  of  the 
observed  mixing  at  the  edge  of  turbulent  combustion  layers.  The 
thickness  of  this  PCL  in  the  model  was  smaller  than  l  mm,  as 
shown  in  Figure  6,  This  figure,  which  shows  the  value  of  the 
dimensionless  gravity  coefficient  g'(z)  in  the  two  circular 
models  also  confirms  the  validity  of  Eq.  2.  Near  the  ground; 
however,  the  value  of  g'  may  be  affected  by  its  initial  value  at 
the  source. 

Since  g'  is  proportional  to  Ap/pa,  which  is  approximately 
proportional  to  AT/Ta,  the  order  of  magnitude  of  AT/Ta  in  large- 
area  fires  was  estimated  (Figure  7). 

The  conclusion  that  AT/T&  in  the  fire  plume  decreases  with 
the  size  of  the  fire  appears  surprising,  but  it  is  a  direct 
consequence  of  Eq.  2.  One  may  also  conclude  that  the  aY££3S£ 
oxygen  depletion  in  large-area  fire  plumes  is  smaller  than  in 
small-area  fires.  This  conclusion  does  not  imply  that  large 
oxygen  depletion  cannot  exist  within  individual  fires  burning 
within  a  LAF. 

induced  Radial  Velocities 

The  induced  radially- inward  velocities  at  the  edge  of  the 
fire  were  found  to  be 

UR  -  k(oR)1/3 

where  k  *  0.24  -  0.5,  much  smaller  than  in  numerical  models  and 
in  reports  on  large-area  fires.  (It  should  be  stressed  that  the 
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above  estimate  is  for  fires  in  a  calm,  dry,  neutrally  stratified 
atmosphere. ) 

Figure  8  compares  the  shape  of  the  radial  velocity  profile 
for  a  point  source  of  buoyancy  (a)  and  for  an  area  source  of 
buoyancy  (b),  which  can  be  described,  at  z  >  R,  as  a  convective 
plume  from  a  virtual  source  at  Zy.  The  figure  demonstrates  why, 
for  the  same  buoyancy  flux,  the  induced  radial  velocities 
increase  when  zv  is  decreased.  The  conclusion  is  that  area 
sources  produce  large  inward-radial  velocities  which  depend  on 
the  effective  necking  of  the  plume. 

Estimated  Plume_Rise_ln  a  Stablv- St ratified  Atmosphere 

Briggs  equation  or  the  analysis  of  Horton,  Turner  and  Taylor 
are  for  a  buoyant  plume  from  a  point  source.  The  initial  mixing 
in  the  radial  convective  boundary  layer  ( CBL )  of  an  area  source 
would  decrease  the  final  plume  rise. 

Figure  9  shows  a  rough  estimate  of  the  plume  rise  for  an 
area  source  with  a  given  total  buoyancy  flux  as  a  function  of  the 
position  of  the  virtual  origin.  The  estimate  is  based  on  the 
assumption  that  the  .behavior  of  the  plume  at  z  >  0  may  be 
calculated  using  the  width  of  the  virtual  plume  at  z  =  0.  The 
data  is  for  a  buoyancy  flux  produced  by  a  100  km2  fire. 

The  average  plume  rise  for  zy  =  o  is  identical  with  the 
estimate  given  in  the  NCR  report  (1985,  p.  74).  One  clearly  sees 
from  Figure  10  that  the  plume  rise  is  reduced  drastically  with 
the  lowering  of  the  virtual  origin.  At  z y/R  *  -1.5  one  gets  an 
average  plume  rise  of  the  order  of  5  km. 
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This  drastic  reduction  is  initially  surprising.  However,  if 
one  examines  the  theoretical  solution  of  Morton  et  al.  for  the 
loss  of  buoyancy  flux  from  a  point  source  with  height  (Figure 
10),  which  shows  that  most  of  the  buoyancy  flux  is  lost  when  the 
plume  achieves  a  large  width,  the  reason  for  the  large  difference 
in  plume  rise  between  point  sources  and  area  sources  becomes 
clear. 

It  is  interesting  to  note  that  while  lowering  of  zv 
decreases  the  plume  rise,  it  increases  the  induced  radial 
velocities. 


The  CBL  is  the  largest  region  above  the  combustion  layer. 
Our  study  indicates  that  it  is  very  similar  to  atmospheric  CBL: 

•  The  Monin-Obukhov  length  is  very  small  L/R  =  odo-4). 

•  The  effect  of  the  horizontal  shear  is  small  and  the 
induced  horizontal  velocities  do  not  cause  mixing. 

•  Buoyancy  flux  is  in  the  form  of  thermals  with  large 
upward  velocities,  w  >  uD. 

•  To  satisfy  continuity  fresh  cool  air  descends  to  the 
ground  between  the  thermals  (or  the  individual  fire 
plumes) . 

•  Rapid  dilution  of  individual  fire  plumes  is  expected. 

The  flow  visualization  in  Figure  ll  shows  the  radial  CBL  in 
the  two-dimensional  fire  model.  Figure  12  shows,  as  well,  a 
smoke  filament  introduced  outside  the  fire. 

Figure  13  shows  smoke  filaments  introduced  from  outside  the 
fire  as  they  enter  the  CBL.  A  rapid  descent  of  the  fresh  air  is 


195 


.  ^  .  V.  A  1.  «*.  «■.  .vV.V  .V 


j-'v.-  j -‘Sv-v-; 


observed,  as  in  cases  of  fumigation.  The  same  pattern  was 
observed  in  the  circular  models  and  in  Project  Flambeau.  Many 
numerical  and  analytical  models  neglect  this  unique  form  of 
mixing. 


The  Effect  of  Urban  Streets 
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the  flow  visualization  to  be  unchanged. 

Conclusions 

•  Physical  models  can  simulate  important  fe_ures  of  large 
area  fire  plumes. 

•  The  rise,  in  a  stably  stratified  atmosphere,  of  a  plume 
with  a  given  buoyancy  flux  decreases  with  the  area  of  the 
source. 

•  Area  sources  of  buoyancy  induce  larger  radial,  inward, 
ground-level  velocities  than  point  sources,  but  the 
magnitude  of  these  velocities  is  not  as  high  as 
previously  reported. 

•  The  interaction  between  simultaneously  burning  fires  over 
a  large  area  is  not  expected  to  be  large.  This 


conclusion  is  based  on  two  observations:  the  average 
oxygen  depletion  in  plumes  of  LAF  is  not  expected  to  be 
large.  The  upward  velocities  of  the  individual  plumes 


within  a  LAF  are  expected  to  be  larger  than  the  induced 
horizontal  velocities. 

Proposed  Future  Simulations 

•  Measure,  in  small-scale  models,  the  instantaneous 
temperature  (oxygen  depletion)  LAF  and  in  an 
individual  fire  surrounded  by  other  fires. 

•  Study  the  detailed  effect  of  streets  and 
noncombustible  areas  on  the  flow  and  temperature 
field  of  a  LAF. 

•  Study  the  behavior  of  LAF  plumes  in  a  stably 
stratified  atmosphere. 

•  Study  the  merger  of  fire  plumes  from  close  area 
fires. 
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Figure  l .  Schematic  description  of  the  experimental 
configurations. 
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Flame  height  correlation  according  to  P.H. 
Thomas  (1963). 
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Figure  8.  Schematic  description  of  induced  radial 
velocities  for  a  point  source  (a)  and  an  area 
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Figure  9.  Estimate  plume  rise  for  a  100  km  fire  in  a 
stably  stratified  atmosphere  as  a  function  of 
the  position  of  the  virtual  source. 
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Figure  10. 


The  relative  reduction  of  buoyancy  flux  F 
with  the  height  z  in  a  stably  stratified 
atmosphere  according  to  Morton  et  al.  (1956). 
(Note  that  most  buoyancy  flux  disappears  at 
large  heights  after  the  plume  has  increased 
its  width. ) 
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Figure  11.  Photographs  of  the  outer  developing  CBL  in 
the  two-dimensional  fire  model. 
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Photographs  of  smoke  filaments  released  above 
the  CBL  in  the  two-dimensional  fire  model. 
(Note  the  rapid  decent  of  the  filaments  in 
the  CBL. ) 
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-  Continuous  thermodynamics 
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EXT  EFFORTS : 


Numerical  Simulation  of 


Lawrence  Livermore  National  Laboratory 


e  numerical  model 


Benchmark 


Height  (km) 


sTTimT 


Sources  of  Differences 


Future  Research 


REGIONAL-SCALE  INTERACTIONS 
BETWEEN  ATMOSPHERIC 
DYNAMICS  AND  AEROSOLS 


Douglas  L.  Westphal 

The  Pennsylvania  State  Univ., 
Department  of  Meteorology, 

Univ.  Park,  PA  16802 

and 

Owen  B.  Toon 

NASA  Ames  Research  Center, 
Space  Sciences  Division, 

Moffett  Field,  CA  94085 


NUCLEAR  WINTER: 
MESOSCALE  TOPICS 


•  Vertical  lofting  of  aerosol  by  the  fire  and  cloud 

•  Coagulation  and  wet  removal  in  the  convec¬ 
tive  plume 

•  Dynamical  perturbation  by  the  fire  and  cloud 

•  Self-generated  bouyancy  of  soot  clouds; 
cloud  formation  in  bouyant  soot  clouds; 
cycling  of  soot  through  cloud  droplets 

•  Spread  of  individual  plumes; 
formation  of  larger,  “GCM-scale”  clouds 

•  Influence  of  diurnal  radiative  forcing  on  the 
response  of  the  atmosphere 

•  Influence  of  geographical  location  and  season 
on  the  above  items 


DAILY- AVERAGED  MINERAL  DUST  CONCENTRATION  (ug/m3) 

AEROSOL  CONCENTRATION  rtf1  ) 


Barbados 
(13°N ,  59°W ) 
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9  JULY,  1974'  TURBIDITY  ANALYSIS  AND  VERTICAL  MOTIONS  FROM  DIABATIC  HEATING  (700mb) 


VERTICAL 

VELOCITY 

SCALE 


1 


I  CMS' 


LONGITUDE  (°W! 
TURBIOITY  CONTOURS 


Turbidity  isopleths  determined  from  satellite  irradiances 
and  diagnosed  700  mb  vertical  motions  associated  with 
radiative  heating  of  dust  for  1200  GMT,  July  9.  1974 
(Karyainpudi,  1979). 


NUMERICAL  MODELS 


Atmospheric  —  Penn  State/NCAR  Limited- 
Area  Model 

•  Primitive  equations,  hydrostatic 

•  Predicted  ground  temperature 

•  Subgrid-scale  vertical  fluxes 

•  Precipitation:  Large-scale  and 
subgrid-scale  (cumulus)  precipitation 

Aerosol  —  Ames  Research  Center’s  Aerosol  Model 

•  Continuity  equation  for  each  size  aerosol  and 
each  type  of  aerosol 

•  Dry  and  wet  removal 

•  Coagulation 

•  Condensation 

•  Surface  or  gridpoint  source,  or  initial  cloud 


•  Radiation 


x">narKxwxKM  nx-jwj iww.’j  VW.7_V.v;y.V.'  VJ  y.'lfj  V*  T* !.'  7*  »7T,*  V  TV."-*?'* 1  pr 1  ,■>  ',’<  y "  ■ "  y"  *rw7wTwT 

SCAVENGING  PARAMETERIZATION 

Scavenging  rate  A  in  (s'-1)  is  given  by 

A  =  4.2  x  10 ~4ER0-79  R0  6Z  <  7000/EH 
A  =  3R016/H  R0-63  >  7000/.Etf 

where: 

E  is  the  collection  efficiency,  0.83 
U  the  depth  of  the  precipitating  cloud  (m) 

R  the  rainfall  rate  ( mmhr~l ) 


Trrrrrrr*?*, 


-« 

J3 

I 

4 


n 


•  Marshall-Palmer  droplet  size  distribution  as¬ 
sumed 

•  No  dependence  on  aerosol  size 

•  Non-convective  rainfall  Rn  is  due  to  super¬ 
saturation  in  the  grid  box 

•  Convective  rainfall  Rc  and  the  fraction  ap  cov¬ 
ered  by  convection  are  determined  by  a  cu¬ 
mulus  cloud  parameterization 

The  total  scavenging  rate  for  a  grid  box  is 
A  =  (1  —  ap)A(Rn)  +  apA(Rn  4-  Rc ) 


A 


a 
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FUTURE  WORK 

Include  radiative  heating  in  the  Saharan  dust 
simulations 

Conduct  tests  to  determine  the  resolution  require 

ments  of  soot  cloud  modeling 

Use  a  non-hydrostatic  model  to  develop  an 
input 

dataset  of  the  initial  spatial  and  size  distri¬ 
bution  of  a  soot  cloud 

•  Carry  out  experiments  that  address  the  out¬ 
standing  problems  of  nuclear  winter 


Overview  of  Ames  Program 


I  Develop  Natural  Analogs  to  Nuclear  Winter 

A.  Saharan  Dust  Storms 
D.  Yolcanic  Clouds 

II  Determine  Fuel  Loading 

III  Develop  Algorithms  for  Aerosol  Physics 

IV  Perform  Regional  Scale  Modeling 
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Algorithms  for  Aerosol  Physics 


I  Aerosol  Microphysics  and  Transport 


II  Radiative  Transfer  in  Aerosols 


III  Interactions  Between  Aerosols  and  Clouds 


Aerosol  Microphysics  and  Transport 

I  General  Properties 

A.  In  Arbitrary  Coordinates  So  Can  Be  Used 
In  A  Variety  of  Dynamical  Models 

B.  Description  in  Preparation 

II  Transport 

A  Timesplit  So  Can  Be  Used  in  1,2, or  3-D  Models 

B.  Finite  Elements  in  Horizontal,  Curve  Fitting 
Scheme  in  Vertical.  Exactly  Conserve  Mass,Have 
Little  Numerical  Diffusion,No  Negative  Numbers, 
Numerically  Efficient 

III  Microphysics 

A.  Resolves  Size  Distribution 

B.  Includes.Sedimentation,  Coagulation,Multiple 
Aerosol  Constituents,Condensational  Growth, 
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TABLE  1 


Conversion  scaling  factors  for  various  coordinate  systems 


.Coordinate 

U 

!».j 

H  „j 

dsl 

ds^ 

''-Systems-'^ 

Spheri c 

s  »n  3? 

1 

a 

Aj 

Coordinates 

Mercator 

cos  $ 
cos  i-6 

_  | 

cos  A 
cos  4. 

1 

l 

dx 

dy 

Polar 

Stereographi c 

H f 9  l 

1  +  sin  f 
2 

-  1 

-9  1 

dx 

dy 

Conformal 

(l  +  sinf)cos<P» 
(l  +  sin</^)cos<p 

(jt-S<o>})  C*><k  1 

UVSVn^CASC, 

dx 

dy 

Rectangular 

1 

1 

1 

dx 

dy 

Pressure 

Sigma 

P-Ptop  =  P-Ptop 
PUm-Ptop  P* 

Log  Pressure  . 
Z  -  H  In  (P/Psfc) 


1/S 

p*/ 


u  =  Uc/m^v  -_v^2  w  =  ds3/dt  +  Vfaii  V(n-1  c  =  Czvm  H^v. 

KX1  =  ]/(i;K  Kx2  =  V&X  ,  KJ(3  =  Vm"2 

=  S  Ws  dxl  =&"gds  i  dx2  =tt«*3  ds2  dx3  =  ds3 

4  =  latitude  ~  latitude  at  which  projection  is  true 

>.  =  longitude  n  =  constant  which  can  be  adjusted  to  make£»_^ 
a  =  radius  of  Farth(£he  projection  tree  at  a  second  1  at i tude  fi") 

^  =  0.\rS  n  -  (k(k-l)’l  k  =1  n  l  cosfol  -1  In  (t.\n  (rt/y  -  P  o/2  *) 

J  fcosif,  '  '  tanfrr  /4  -  <£,/2  } 

g  =  gravitational  <xCc  <*Wwt  m  a  H  =  Sa»t  height  =  RTo/g 
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AD-A18S  15# 
UNCLASSIFIED 


TECHNICAL  PAPERS  PRESENTED  AT  THE  DEFENSE  NUCLEAR 
AOENCV  GLOBAL  EFFECTS  R.  .  <U>  DOD  NUCLEAR  INFORHATION 
AND  ANALVSIS  CENTER  SANTA  BARBARA  CA. .  15  HAV  BE 

DASIAC-TN-8C-29-V0L-2  DNA881-82-C-8274  F/Q  15/C.  4 


GEOMETRICALLY  INCREASING  VOLUME  BINS: 
v^-'v.Av^-'vf  V, 


-j  K- REFERENCE  BIN  WIDTH,  V0«y  V, 


0  12  3  4  5 

DEAD  „  BIN 

VOLUME  u 

'CLASSICAL'  VOLUME  BIN  SEQUENCE: 

Vj « i  V,;  AVj*V0«V, 

BIN  CENTRAL  BIN 
-j  J— BIN  WIDTH,  V0  BOUNDARIES  VOLUME,  Vj 

■A  i  i  i  i  i  i  i  i  i  i  i  r'T'i  i  i  i  i  \  i  i 

a  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i  i\  i  i 
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Figure  7.-  Aerosol  particle  volume  bin  structures  for  numerical  models. 
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Figure  8.-  Aerosol  droplet  growth.  The  evolution  of  a  droplet  size 
distribution  following  a  change  in  the  droplet  growth  rate  is  shown 
Both  analytical  and  numerical  solutions  are  given.  The  conditions 
for  the  calculation  are  described  in  the  text. 
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Figure  9.-  Aerosol  coagulation.  Calculated  transient  decay  curves  are 
shown  for  the  number  of  particles  in  the  smallest  model  size  bin, 
and  for  the  total  number  of  particles  of  all  sizes,  when  one  particle 
is  initially  placed  in  the  smallest  model  bin.  Time  is  measured 
relative  to  the  characteristic  coagulation  time,  t  .  The  exact 
solutions  of  Smoluchowski  are  shown  for  comparison. 

PARTICLE  RADIUS,  ftm 
0.01  Ql  1.0  3.0 


Figure  10.-  The  development  of  an  aerosol  size  spectrum  by  coagulation. 
The  conditions  for  the  calculation  are  the  same  as  for  figure  9. 
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Radiative  Transfer 


I  Solar  Wavelengths 

Use  Two-Stream  Approximation  for  Scattering  with 
Tridiagonal  Solver.  Numerically  Efficient,  Accurate 
to  10%  or  Better  for  Heating  Rates. 


II  Infrared  Wavelengths 

Use  Two  Stream  Source  Function  Technique  for 
Scattering.  Numerically  Efficient,  Exact  in  Limit  of 
No  Scattering,  10%  Accuracy  with  Scattering. 

III  Status 

A.  Algorithms  Developed  and  Checked  for  Speed  and 
Accuracy 

B.  Absorption  Coefficients  Under  Development 
Based  on  McClatchey  tape. 
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F+  =  klnexp{\nr)  +  Tnkinexp{-\nr)  +  <7+  (r)  (18) 

*n  =r„*in«*p(Anr)  +  *jne*p(-A,r)+C"(r).  (19) 
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Mesoscale  Modeling  of  Coastal  Flows 


Numerical  Simulation  Using  the  Colorado  State  University 
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CSU  Mesoscale  Model  —  Radiation  Calculation 
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MODIFICATIONS  TO  THE  CSU  MESOSCALE  MODEL 
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INHERENT  ASSUMPTIONS  FOR  NUCLEATION  SCAVENGING 
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ATTN:  AWARSHAWSKY 
ATTN.  J  WANG 
ATTN:  WSHIH 

PHYSICAL  RESEARCH  INC 
ATTN:  R  JORDANO 

PHYSICAL  RESEARCH,  INC 
ATTN:  DWESTPHAL 
ATTN:  DWHITENER 
ATTN.  H  WHEELER 
ATTN  R  BUFF 
ATTN:  R  DELIBERIS 
ATTN:  T  STEPHENS 
ATTN:  WC  BLACKWELL 

PHYSICAL  RESEARCH.  INC 
ATTN:  G HARNEY 
ATTN:  J  DEVORE 
ATTN:  J  THOMPSON 
ATTN:  RSTOECKLY 
ATTN  W  SCHLEUTER 

PHYSICAL  RESEARCH,  INC 
ATTN:  H  SUGIUCHI 

POLYTECHNIC  OF  NEW  YORK 
ATTN.  BJ  BULKIN 
ATTN:  GTESORO 

PRINCETON  UNIVERSITY 

ATTN:  J  MAHLMAN 

QUADRI  CORP 

ATTN:  HBURNSWORTH 


R  &  D  ASSOCIATES 

ATTN:  AKUHL 
ATTN:  F  GILMORE 
ATTN:  G  JONES 
ATTN:  J SANBORN 
ATTN:  RTURCO 

R  4  0  ASSOCIATES 

ATTN:  B  YOON 

R  J  EDWARDS  INC 

ATTN:  R  SEITZ 

RADIATION  RESEARCH  ASSOCIATES,  INC 
ATTN:  B  CAMPBELL 
ATTN:  M  WELLS 

RAND  CORP 

ATTN.  GL  DONOHUE 
ATTN:  P  ROMERO 

RAND  CORP 

ATTN:  J  GERTLER 

ROCKWEL.  INTERNATIONAL  CORP 
ATTN:  SI  MARCUS 

ROCKWELL  INTERNATIONAL  CORP 
ATTN:  J  KELLEY 

S-CUBED 

ATTN:  B  FREEMAN 
ATTN:  KDPYATT,  JR 
ATTN:  RLAFRENZ 

S-CUBED 

ATTN:  C NEEDHAM 
ATTN:  SHIKIDA 
ATTN:  T  CARNEY 

SCIENCE  APPLICATIONS  INC 
ATTN:  REDELMAN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  CHILL 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  D  HAMLIN 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  B  MORTON 
ATTN:  B  SCOTT 
ATTN:  D  SACHS 
ATTN:  GT  PHILLIPS 
ATTN:  JBENGSTOM 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  D  BACON 
ATTN:  DRLGOURE 
ATTN:  FGIESSLER 
ATTN:  J COCKAYNE 
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ATTN:  J  SHANNON 
ATTN:  J STUART 
ATTN:  MSHARFF 
ATTN:  WLAYSON 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN.  JSONTOWSKI 

SCIENCE  APPLICATIONS  INTL  CORP 
ATTN:  T  HARRIS 

SCIENTIFIC  RESEARCH  ASSOC,  INC 
ATTN:  B  WEINBERG 

SPARTA  INC 

ATTN:  R HARPER 

SRI  INTERNATIONAL 

ATTN.  CWITHAM 
ATTN:  O  GOLDEN 
ATTN:  D  MACDONALD 
ATTN.  D  ROBERTS 
ATTN:  EUTHE 
ATTN:  GABRAHAMSON 
ATTN:  J  BACKOVSKY 
ATTN:  W  CHESNUT 
ATTN:  W  JOHNSON 

SRI  INTERNATIONAL 

ATTN.  RBRAMHALL 
ATTN:  RWOOLFOLK 
ATTN:  W  VAIL 

STAN  MARTIN  ASSOCIATES 
ATTN:  SB  MARTIN 

STANTON  CONSULTING 

ATTN:  M STANTON 

SWETL,  INC 

ATTN:  TY  PALMER 

SYSTEM  PLANNING  CORP 
ATTN.  J  SCOURAS 
ATTN:  MBIENVENU 
ATTN.  R  SCHEERBAUM 

SYSTEMS  AND  APPLIED  SCIENCES  CORP 
ATTN:  M  KAPLAN 

TECHNOLOGY  INTERNATIONAL  CORP 
ATTN:  WBOQUIST 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  0  ORMOND 
ATTN:  F LEOPARD 
ATTN.  J  FORD 

TELEDYNE  BROWN  ENGINEERING 
ATTN:  DGUICE 
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TEXAS  ENGR  EXPERIMENT  STATION 
ATTN:  WH  MARLOW 

TOYON  RESEARCH  CORP 
ATTN.  CTRUAX 
ATTN:  J  GARBARINO 
ATTN:  JISE 


TRW 

ATTN.  HBURNSWORTH 
ATTN:  JBELING 

TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN.  FFENDELL 
ATTN:  G  KIRCHNER 
ATTN:  GMR02 
ATTN:  H  CROWDER 
ATTN:  JFEDELE 
ATTN:  M  BRONSTEIN 
ATTN.  R  BACHARACH 
ATTN:  S  FINK 
ATTN  T  NGUYEN 


TRW  ELECTRONICS  &  DEFENSE  SECTOR 
ATTN:  M  HAAS 


VISIDYNE,  INC 

ATTN.  H  SMITH 
ATTN:  J CARPENTER 

WASHINGTON,  UNIVERSITY  OF 
ATTN:  J  I  KATZ 

FOREIGN 

AERE  ENVIRONMENTAL  AND  MEDICAL  SC 
ATTN  SPENKETT 

ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT 
ATTN:  PFARICHARDS 

ATOMIC  WEAPONS  RESEARCH  ESTABLISHMENT 
ATTN:  D  L  JONES 
ATTN  DM  MOODY 

AUSTRALIA  EMBASSY 

ATTN.  DR  LOUGH 

ATTN:  MAJ  GEN  H  J  COATES 

ATTN:  P  PROSSER 

BRITISH  DEFENCE  STAFF 
ATTN:  C  FENWICK 
ATTN:  J  CRANIDGE 
ATTN:  J  EDMONDS 
ATTN  M  NORTON 
ATTN:  P  WEST 

CANADIAN  FORESTRY  SERVICE 
ATTN:  B  STOCKS 
ATTN:  TLYNHAM 
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ATTN:  IGALBALLY 


CSIRO:  ATMOSPHERIC  RESEARCH 
ATTN:  A  PITTOCK 


EMBASSY  OF  BELGIUM 

ATTN:  LARNOULO 


ISRAEL  EMBASSY 

ATTN:  N  BELKIND 


MAX  PLANCK  INSTITUTE  FOR  CHEMISTRY 
ATTN:  PJCRUT2EN 


MINISTRY  OF  DEFENCE 
ATTN:  R  RIDLEY 


NATIONAL  DEFENCE  HEADQUARTERS 
ATTN:  HAROBITALLE 


TRINITY  COLLEGE 

ATTN:  F  HARE 


DIRECTORY  OF  OTHER 


ATMOS.  SCIENCES 

ATTN:  GSISCOE 


BROWN  UNIVERSITY 

ATTN:  RK  MATTHEWS 


BUCKNELL  UNIVERSITY 

ATTN:  0  ANDERSON 


CALIFORNIA,  UNIVERSITY 

ATTN:  R  WILLIAMSON 


CALIFORNIA,  UNIVERSITY  OF 

ATTN:  L  BADASH/DEPT  OF  HISTORY 


COLORADO,  UNIVERSITY  LIBRARIES 
ATTN:  JBIRKS 
ATTN:  R  SCHNELL 


DREXEL  UNUVERSITY 

ATTN:  J  FRIEND 


DUKE  UNIVERSITY 

ATTN:  F  DELUCIA 


GEORGE  MASON  UNIVERSITY 
ATTN:  PROFS  SINGER 
ATTN:  R  EHRLICH 


GEORGE  WASHINGTON  UNIVERSITY 
ATTN:  R  GOULARD 


GEORGIA  INST  OF  TECH 

ATTN:  E  PATTERSON 
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HARVARD  COLLEGE  LIBRARY 
ATTN:  W  PRESS 


HARVARD  UNIVERSITY 

ATTN:  G  CARRIER 


HARVARD  UNIVERSITY 

ATTN:  DEARDLEY 


IOWA,  UNIVERSITY  OF 

ATTN:  HISTORY  DEPT/S  PYNE 


MARYLAND  UNIVERSITY  OF 

ATTN:  A  ROBOCK  DEPT  METEOROLOGY 
ATTN:  A  VOGELMANN  DEPT  METEOROLOGY 
ATTN:  R  ELLINGSON  DEPT  METEOROLOGY 


MIAMI  LIBRARY  UNIVERSITY  OF 
ATTN:  C  CONVEY 


MIAMI  UNIV  LIBRARY 

ATTN:  J  PROSPERO  ATMOS  SC 


NEW  YORK  STATE  UNIVERSITY  OF 
ATTN:  R  CESS 


OAK  RIDGE  ASSOCIATED  UNIVERSITIES 
ATTN:  C  WHITTLE 


PENNSYLVANIA  STATE  UNIVERSITY 
ATTN:  DWESTPHAL 


SOUTH  DAKOTA  SCH  OF  MINES  *.  TECH  LIB 
ATTN:  H  ORVILLE 


TENNESSEE,  UNIVERSITY  OF 
ATTN.  K  FOX 


UNIVERSITY  OF  SOUTH  FLORIDA 
ATTN:  SYING 


UNIVERSITY  OF  WASHINGTON 


CLEOVY 
L  RAOKE 
P  HOBBS 


VIRGINIA  POLYTECHNIC  INST  LIB 
ATTN:  MNADLER 


WASHINGTON  STATE  UNIVERSITY 
ATTN:  DR  A  CLARK 


WISCONSIN  UNIVERSITY  OF 
ATTN:  P  WANG 
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